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Abstract 

The emergence of rotational bands is observed in no-core configuration interaction (NCCI) calculations for the odd-mass Be iso- 
(f) topes (7 < A < 13) with the JISP16 nucleon-nucleon interaction, as evidenced by rotational patterns for excitation energies, 
<~H quadrupole moments, and E2 transitions. Yrast and low-lying excited bands are found. The results demonstrate the possibility of 

Well-developed rotational structure in NCCI calculations using a realistic nucleon-nucleon interaction. 
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1. Introduction 

Nuclei exhibit a wealth of collective phenomena, including 
clustering, rotation, and pairing [1-3]. Collective dynamics 
have been extensively modeled in phenomenological descrip- 
tions [1, 3-5]. Some forms of collectivity may also be obtained 
J microscopically in the conventional (valence) shell model, e.g., 
Elliott SU(3) rotation [6, 7]. However, observing the emergence 
of collective phenomena directly from first principles — that is, 
in a fully ab initio calculation of the nucleus, as a many-body 
system in which all the constituent protons and neutrons partic- 
ipate, with realistic interactions — remains as an outstanding 
challenge. 

Recent developments in large-scale calculations have 
brought significant progress in the ab initio description of 
light nuclei [8-12]. In ab initio no-core configuration in- 
teraction (NCCI) approaches — such as the no-core shell 
model (NCSM) [11, 13-16], no-core Monte Carlo shell model 
(MCSM) [17], or no-core full configuration (NCFC) [18] meth- 
ods — the nuclear many -body bound-state eigenproblem is for- 
mulated as a matrix diagonalization problem. The Hamilto- 
nian is represented with respect to a basis of antisymmetrized 
products of single-particle states, generally harmonic oscillator 
states, and the problem is solved for the full system of A nu- 
cleons, i.e., with no inert core. In practice, such calculations 
must be carried out in a finite space, obtained by truncating 
the many-body basis according to a maximum allowed number 
N max of oscillator excitations above the lowest oscillator con- 
figuration (e.g., Ref. [11]). With increasing N mlLX , the results 
converge towards those which would be achieved in the full, 
infinite-dimensional space for the many-body system. 

Computational restrictions limit the extent to which con- 
verged calculations can be obtained for the observables needed 
for the identification of collective phenomena. In particular, 
the observables most indicative of rotational collectivity — E2 



matrix elements — present special challenges for convergence 
in an NCCI approach [19, 20], due to their sensitivity to the 
large-radius asymptotic portions of the nuclear wave function. 
Nonetheless, some promising suggestions of collective phe- 
nomena, e.g., deformation and clustering, have already been 
obtained in ab initio calculations [20-25]. 

In this letter, we observe the emergence of collective rotation 
in ab initio NCCI calculations for the Be isotopes, using the 
realistic JISP16 nucleon-nucleon interaction [26]. Evidence for 
rotational band structure is found in the calculated excitation 
energies, quadrupole moments, and E2 transition matrix ele- 
ments. In NCCI calculations of the even-mass Be nuclei, yrast 
sequences of angular momenta 0, 2, 4, . . . arise with calculated 
properties resembling those of K = ground-state rotational 
bands (see Ref. [27] for a preliminary report of comparable re- 
sults for 12 C). However, the most distinctive, well-developed, 
and systematic rotational band structures are observed in calcu- 
lations for odd-mass nuclei. Given the same range of excitation 
energies and angular momenta, the low-lying AJ = 1 bands in 
the odd-mass nuclei provide a richer set of energy and electro- 
magnetic observables. We therefore focus here on the odd-mass 
Be isotopes, specifically, with 7 < A < 13. After a brief re- 
view of the properties expected in nuclear rotational structure 
(Sec. 2), the results for rotational bands in NCCI calculations 
of these Be isotopes are presented (Sec. 3). Preliminary results 
for 9 Be were reported in Ref. [25]. 



2. Rotation 

We first review nuclear collective rotation and its expected 
signatures [1, 3]. Under the assumption of adiabatic separa- 
tion of the rotational degree of freedom, a nuclear state may be 
described in terms of an intrinsic state, as viewed in the non- 
inertial intrinsic frame, together with the rotational motion of 
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this intrinsic frame. For axially symmetric structure, in partic- 
ular, the intrinsic state \<p K ) is characterized by definite angular 
momentum projection K along the intrinsic symmetry axis. The 
full nuclear state \\Itjkm), with total angular momentum J and 
projection M, then has the form 

+{-) J+K 9 j M _ K {n<Pk\m (i) 

where § represents the Euler angles for rotation of the intrinsic 
state, and is the ^-conjugate intrinsic state, which has 
angular momentum projection -K along the symmetry axis. 

The most recognizable features in the spectroscopy of rota- 
tional states reside not in the states taken individually but in 
the relationships among the different states \i)/jkm) sharing the 
same intrinsic state \(f> K ). These states constitute members of 
a rotational band, with angular momenta J = K, K + 1, . . ., 
except with only even J (or only odd J, depending upon the 
intrinsic &2 symmetry) for K = bands. Within a rotational 
band, energies follow the pattern 

E(J) = E Q +AJ(J + 1), (2) 

where, in terms of the moment of inertia J about an axis per- 
pendicular to the symmetry axis, A = h 2 /(2ST). For K = 1/2 
bands, the Coriolis contribution to the kinetic energy results in 
an energy staggering given by 

E(J) = E + A[J(J + 1) + a{-) J+ll \j + \)\ (3) 

where a is the Coriolis decoupling parameter. Reduced ma- 
trix elements {^j f K\\Q2\\^j,K) of the electric quadrupole oper- 
ator Q2 between states within a band are entirely determined 
by the rotational structure, except for the overall normaliza- 
tion, which is proportional to the intrinsic quadrupole moment 
eQo = (167t/5) 1/2 ((/>a:|22,oI<Aa:)- In particular, quadrupole mo- 
ments within a band are obtained as 

and reduced transition probabilities as 

B(E2; Ji -» J f ) = -^{JiKltyJfKfieQof. (5) 
ion 

In obtaining these results, Q2 can be taken to be any operator of 
the form Q2 M = Y*t=i e irf ;) and may therefore represent the 
electromagnetic E2 transition operator, mass quadrupole tensor, 
neutron quadrupole tensor, etc., depending upon the choice of 
coefficients e, (see Sec. 3). 

3. Results 

An NCCI calculation is defined by the interaction for the nu- 
clear system and by the truncated many-body space in which 
the calculation is carried out. The present calculations use the 



JISP16 interaction [26], which is a two-body interaction de- 
rived from neutron-proton scattering data and adjusted via a 
phase-shift equivalent transformation to describe light nuclei 
without explicit three-body interactions. The bare interaction 
is used, without renormalization to the truncated space [18]. 
The Coulomb interaction has been omitted from the Hamilto- 
nian, to ensure exact conservation of isospin, thereby simplify- 
ing the spectrum. (The primary effect of the Coulomb interac- 
tion, when included, is to induce a shift in the overall binding 
energies, which is irrelevant to identification of rotational band 
structure.) These calculations are carried out for oscillator trun- 
cations ranging from N m - dx = 10 for 7 Be to N max = 7 for 13 Be, 
with basis oscillator HQ. parameters near the variational min- 
imum (HQ. = 20-22.5 MeV). The proton-neutron M-scheme 
code MFDn [28-30] has been used for the many-body calcula- 
tions. 

The calculated excitation energies for low-lying states of the 
odd-mass Be isotopes (with 7 < A < 13) are shown in Figs. 1 
and 2. For each nucleus, there are two parity spaces to consider, 
shown separately in these two figures (energies are taken rela- 
tive to the lowest state of the same parity). We refer to the par- 
ity of the lowest allowed oscillator configuration (negative for 
7 ' 9 ' U Be, positive for 13 Be) as the natural parity (Fig. 1) and that 
obtained by promoting one nucleon by one shell as the unnat- 
ural parity (Fig. 2). The NCCI bases for these spaces consist 
of states with even and odd numbers of oscillator excitations, 
respectively, above the lowest configuration. While the low- 
est unnatural parity states normally lie at significantly higher 
energy than those of natural parity, they are calculated to lie 
within a few MeV of the lowest natural parity states in the iso- 
topes 9 ' 1113 Be [25], which are therefore included in Fig. 2. Note 
that parity inversion arises for 11 Be, i.e., the ground state is ex- 
perimentally [31] in the unnatural parity space, and both spaces 
are near-degenerate in calculations at finite N m - dx (see Ref. [11]). 
The minimal isospin (T = T z ) spectrum is shown in each case. 

To facilitate identification of rotational bands, it is helpful to 
plot the calculated excitation energies with respect to J( J +1), 
so that energies within an ideal rotational band would lie on 
a straight line — or, for K = 1/2 bands, staggered about a 
straight line. For the candidate K = 1/2 bands in Figs. 1 and 2, 
an energy fit is obtained by adjusting the parameters of (3) to 
the first three bandmembers (the remainder of the line is thus an 
extrapolation). For the remaining bands, a straight line fit of (2) 
to all bandmembers is shown. 

The yrast and near-yrast states yield the most immediately 
recognizable sets of candidate bandmembers. Yrast rotational 
bands (with bandmembers indicated by solid black squares in 
Figs. 1 and 2) are found with K = 1/2 except in the natural 
parity space of 9 Be [Fig. 1(b)], for which the yrast band has 
K = 3/2. The density of states rapidly increases off the yrast 
line, hindering identification of candidate bands and further- 
more suggesting that the rotational states may be fragmented 
by mixing with nearby states. Nonetheless, several excited can- 
didate bands (indicated by solid red squares) can also be clearly 
identified, with 1/2 < K < 5/2, once transition strengths have 
been taken into account. 

For the yrast K = 1/2 bands, as a result of Coriolis decou- 
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Figure 1: Excitation energies obtained for states in the natural parity spaces 
of the odd-mass Be isotopes: (a) 7 Be, (b) 'Be, (c) "Be, and (d) 13 Be. Ener- 
gies are plotted with respect to J(J + 1) to facilitate identification of rotational 
energy patterns, while the J values themselves are indicated at top. Filled sym- 
bols indicate candidate rotational bandmembers (black for yrast states and red 
for excited states). The lines indicate the corresponding best fits for rotational 
energies. Where quadrupole transition strengths indicate significant two-state 
mixing (see text), more than one state of a given J is indicated as a bandmem- 
ber. 

pling, it should be noted that alternate bandmembers are raised 
in energy into a region of higher density of states, which com- 
plicates identification and is conducive to fragmentation. The 
energy staggering in the calculated yrast band of the 7 Be nat- 
ural parity space [Fig. 1(a)] — in which the 7 = 3/2, 7/2, 
... levels are lowered, and the 7 = 1/2, 5/2, ... levels are 
raised — corresponds to a negative value of the decoupling 
parameter. Note that the staggering is sufficiently pronounced 
that the two lowest-7 bandmembers are inverted, as experimen- 
tally observed for this nucleus [32]. Then, positive values of 
the decoupling parameter are instead obtained for the remain- 
ing K - 1/2 bands. 

It is interesting to compare these ab initio results for the yrast 
bands in the natural parity spaces (Fig. 1) with the Nilsson 
model predictions [33, 34]. Specifically, the calculated yrast 
bands have K = 1/2 {a « -1.4) for 7 Be, K = 3/2 for 9 Be, 
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Figure 2: Excitation energies obtained for states in the unnatural parity spaces 
of the odd-mass Be isotopes: (a) 'Be, (b) "Be, and (c) I3 Be. Energies are 
plotted with respect to J(J + 1 ) to facilitate identification of rotational energy 
patterns, while the J values themselves are indicated at top. Filled symbols 
indicate candidate rotational bandmembers (black for yrast states and red for 
excited states). The lines indicate the corresponding best fits for rotational 
energies. Quadrupole transition strengths indicate significant but ambiguous 
two-state or multistate mixing for certain poorly-isolated bandmembers. 



K= 1/2 (a * +1.2)for n Be, andtf= 1/2 (a « +3.1) for 13 Be. 
The expected Nilsson [Nn z A£l] asymptotic quantum number 
assignments (see Fig. 5-1 of Ref. [3]) and corresponding Nils- 
son values of the decoupling parameter are [llO^] (a « -1) 
for 7 Be, [101|] for 9 Be, [101±] (a * 0) for n Be, and [220±] 
(a w +1) for 13 Be. We see consistency not only in the K (= Q) 
quantum numbers for the band but also in the qualitative trend 
of the decoupling parameters for these bands. (The Nilsson val- 
ues for a [33, 34] consider mixing of spherical orbitals only 
within a single spherical oscillator shell, which is sufficient for 
a weakly-deformed oscillator-like mean field. However, they 
should not be expected to provide quantitative accuracy for a 
nucleon in, say, the mean field produced by a double-or 8 Be 
core.) 

The quadrupole moments for all states within the candidate 
bands are shown in Figs. 3 and 4, both for the yrast bands 
(black squares) and for the excited bands (red diamonds). The 
values are normalized to Qq, to facilitate comparison with the 
rotational predictions for Q( J)/Qq from (4) (shown as curves 
in each plot). The value of Qq used for normalization has in 
each case been obtained simply from the quadrupole moment 
of the lowest-energy bandmember of nonvanishing quadrupole 
moment. (Thus, for K — 1/2 bands, since the quadrupole mo- 



3 




-0.4 - 7 

(a) Be 




1/2 



3/2 



5/2 



7/2 



9/2 



Figure 3: Quadrapole moments calculated for candidate bandmembers in 
the natural parity spaces of the odd-mass Be isotopes: (a) 7 Be, (b) 9 Be, 
(c) "Be, and (d) 13 Be. The states are as identified in Fig. 1 and are shown 
as black squares for yrast states or red diamonds for excited states. Filled sym- 
bols indicate proton quadrapole moments, and open symbols indicate neutron 
quadrapole moments. The curves indicate the theoretical values for a K = 1/2 
or K = 3/2 rotational band, as appropriate, given by (4). Quadrapole moments 
are normalized to Qg, which is defined by either the J = 3/2 or J = 5/2 band- 
member (see text). 



ment of the 7=1/2 bandhead vanishes identically, either the 
J — 3/2 or 5/2 bandmember is used for normalization, accord- 
ing to the staggering.) Quadrupole moments in Figs. 3 and 4 
are calculated using both the proton (filled symbols) and neu- 
tron (open symbols) quadrupole tensors. 1 

Finally, in-band transition strengths are shown in Fig. 5, 
again as obtained for both proton (solid symbols) and neutron 
(open symbols) quadrupole operators, and for AJ = 2 tran- 
sitions (upper curves, solid) and AJ = 1 transitions (lower 



'The proton quadrapole tensor, defined as Qi^ p = 2jJ=i r 2 ^-^(t p,i), is 
the operator used in calculation of the physically observable electromagnetic 
moments and transitions. However, the rotational relations (4) and (5) are 
equally applicable to matrix elements of the neutron quadrapole tensor, Q2^,n = 
Zjli >i -YtytXnj)- These therefore provide a valuable complementary set of ob- 
servables for purposes of investigating whether or not the nuclear wave func- 
tions satisfy the conditions of adiabatic rotational separation, particularly rel- 
evant, due to the high neutron-proton asymmetry, in the neutron-rich Be iso- 
topes. 
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Figure 4: Quadrapole moments calculated for candidate bandmembers in the 
unnatural parity spaces of the odd-mass Be isotopes: (a) 9 Be, (b) "Be, and 
(c) 13 Be. The states are as identified in Fig. 2 and are shown as black squares for 
yrast states or red diamonds for excited states. Filled symbols indicate proton 
quadrapole moments, and open symbols indicate neutron quadrupole moments. 
The curves indicate the theoretical values for a K = 1/2 or K = 3/2 rotational 
band, as appropriate, given by (4). Quadrapole moments are normalized to Qq, 
which is defined by either the J = 3/2 or J = 5/2 bandmember (see text). 



curves, dashed). The various K - 1/2 bands are superposed in 
Fig. 5(a), and the K = 3/2 bands are shown in Fig. 5(b). Transi- 
tion strengths are normalized as B(E2; J — > J - AJ)/(eQo) 2 , for 
comparison with the rotational values from (5). The same 2o 
values are used as in Figs. 3 and 4, i.e., obtained from 2(3/2) 
or 2(5/2). Therefore, no free normalization parameters re- 
main for the B(E2) strengths in Fig. 5. For instance, it may 
be observed that the values for B(E2;3/2 — > l/2)/(e2o) 2 in 
Fig. 5(a) cluster at the rotational value, indicating that the cal- 
culated B(E2; 3/2 — » 1/2) strengths are in the proper relation 
to the calculated 2(3/2) moment or 2(5/2) moment, as appro- 
priate, consistent with adiabatic rotation. 

The level of resemblance between the calculated energies, 
quadrupole moments, and transition strengths for the candidate 
bands and the expected rotational values in Figs. 1-5, while 
clearly not perfect, indicates a remarkably clean separation of 
rotational and intrinsic degrees of freedom in these ab initio 
NCCI calculations. One should bear in mind that quadrupole 
moments of arbitrarily chosen states in the spectrum fluctuate 
not only in magnitude but also in sign, and that calculated E2 
strengths among arbitrarily chosen pairs of states fluctuate by 
many orders of magnitude. (The 3/2 — » 1/2 transitions in Fig. 5 
are enhanced by factors of ~ 1 . 1-17 relative to the typical Weis- 
skopf single-particle estimate [35].) 

It is worth highlighting a few notable features from the band 
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Figure 5: Transition B(E2; J —> 7—2) and B(E2; J — » J— 1) strengths calculated 
between candidate bandmembers, for bands with (a) K = 1/2 or (b) K = 3/2, 
in the odd-mass Be isotopes 7 Be, 'Be, "Be, and 13 Be (asterisks in the legend 
indicate excited bands). The curves indicate the theoretical values for AJ = 2 
(upper curves, solid) and AJ = 1 (lower curves, dashed) transitions in a band 
of the given K, given by (5). Calculated transition strengths of the same AJ 
within a band are connected by dotted lines. Filled symbols indicate proton 
strengths, and open symbols indicate neutron strengths. Values are normalized 
to the same intrinsic quadrupole moment Qq as in Fig. 3 or 4. Where transition 
strengths indicate significant fragmentation of the rotational state over multi- 
ple calculated states, namely, for the "Be excited J = 7/2 state identified in 
Fig. 1(c), the B(E2) values for transitions into or out of this state are summed. 



structures in Figs. 1-5: 

(1) The K — 1/2 yrast bands in the unnatural parity spaces 
(Fig. 2) can be traced to J values as high as ~ 13/2. For in- 
stance, for 13 Be [Fig. 2(c)], the energies of the J - 7/2, 9/2, 
1 1 /2, and 13/2 bandmembers all agree with the rotational val- 
ues, from (3), to within 0.4 Me V, and a / = 15/2 bandmember 
can also be reasonably identified (within 1 .5 Me V of the ro- 
tational energy). The quadrupole moments [Fig. 4(c)] of the 
J = 3/2 and / = 5/2 bandmembers (the latter is used to de- 
termine Qo in the figure) are in the expected rotational ratio, 
from (4), to within 1.1% for protons or 0.4% for neutrons. The 
quadrupole moments for the higher bandmembers are highly 
consistent between protons and neutrons and have the expected 
sign, but they gradually fall off from the rotational values, ap- 
proaching zero for the J — 15/2 bandmember. 

(2) For the yrast and low-lying rotational bands in the natural 
parity spaces (Fig. 1), rotational behavior appears to terminate 
at generally lower angular momentum. For instance, for 11 Be 
[Fig. 1(c)], the K — 1/2 yrast band terminates at J = 7/2 on 
the basis of energies: the lowest calculated 7 = 7/2 state lies 
within 0.5 MeV of the expected energy extrapolated for an yrast 



bandmember, but the lowest calculated J -9/2 state is 1 1 MeV 
too high in energy to be an yrast bandmember. The terminat- 
ing angular momentum expected in a simple valence /j-shell or 
NCCI N m - dx = description is, in fact, J = 7/2. The quadrupole 
moments [Fig. 3(c)] suggest that the viability of a rotational de- 
scription may end even earlier, at J — 5/2. Similar comments 
may be made about the yrast and excited bands in 7 Be and 9 Be 
[Fig. l(a,b)], where the quadrupole moments [Fig. 3(a,b)] are 
in close agreement with rotational values through 7 = 7/2, but 
then begin to deviate significantly at J = 9/2. 

(3) To some extent in the quadrupole moments, but espe- 
cially in the AJ — 2 transition strengths for the K = 1/2 
bands [Fig. 5(a)], one may observe that the E2 matrix ele- 
ment strengths start at the expected rotational values for low 
/ but then systematically fall off below the rotational values 
at higher /. This trend signals deviation from a strict adia- 
batic rotational picture, as described in Sec. 2, but it is also, at 
least qualitatively, in agreement with more microscopic treat- 
ments of nuclear rotation. Specifically, E2 matrix elements 
within an Elliott SU(3) band decline in strength as band ter- 
mination is approached (see discussion in Ref. [7]). A similar 
falloff can be obtained in Sp(3, R) symplectic calculations [36] 
of rotational bands (see Fig. 6 of Ref. [37]). Whether or not 
such SU(3) or Sp(3, R) rotational mechanisms are at play in the 
present NCCI results awaits full analysis in an SU(3)/Sp(3,R) 
symmetry-adapted implementation of the NCCI approach [38]. 

As we explore the interpretation of NCCI results in a rota- 
tional context, it is interesting to note that straightforward frag- 
mentation of the rotational strength over two calculated levels 
can be observed for the 7 = 7/2 member of the excited band in 
the natural parity space of 1 'Be [Fig. 1(c)]. Transitions into and 
out of this state are fragmented in the approximate proportion 
0.4 : 0.6. However, the summed strengths, shown in Fig. 5(b), 
which combine the fragmented transitions involving this level, 
are in near-perfect agreement with rotational values. 

We also note that a staggering may be observed in the AJ — 2 
transition strengths, to a greater or lesser degree, for the various 
K — 1/2 bands [Fig. 5(a)]. Such staggering is in fact consis- 
tent with the adiabatic rotational picture, once the (ipKlQi.iKl^k) 
cross term in the rotational E2 matrix element, neglected in (5), 
is taken into account [see (6.38) of Ref. [1]]. In well-deformed 
rotor nuclei, this contribution is commonly ignored, on the pre- 
sumption that (<Pk\Q2,o\<Pk) ~ Qo is strongly enhanced while 
(<Pk\Q2,2k\</>k) i s °f single-particle strength [1]. However, in 
light nuclei, where the collective enhancement is weaker, such 
a single-particle contribution may be expected to be nonnegli- 
gible in comparison, and to be of approximately the magnitude 
seen in Fig. 5(a). 



4. Conclusion 

The principal challenge in identifying collective structure in 
NCCI calculations with realistic interactions lies in the weak 
convergence of the relevant observables. Eigenvalues and other 
calculated observables are dependent upon both the truncation 
Nmm and the oscillator length parameter (or HQ.) for the NCCI 
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Figure 6: Relative quadrapole moments Q(J)/Qo, calculated with the (a) proton 
and (b) neutron quadrupole tensor, for the K = 3/2 yrast band in the 'Be natural 
parity space, from calculations with basis truncations N mliK = 4, 6, 8, and 10. 
Since the normalization go is obtained in each case from Q(3/2), the plotted 
values indicate the /V max dependence of the ratio 2(7)/ 2(3/2). 



from SU(3) rotation in an extended multi-shell valence space — 
and the relevance of some form of Nilsson-like strong-coupling 
picture for rotation in ab initio calculations of odd-mass light 
nuclei. Indeed, the proton and neutron density distributions 
found in Ref. [25] for the ground state of 9 Be suggest the emer- 
gence of two alpha clusters, with the additional neutron in a n 
orbital. 
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basis. Although it is possible to extrapolate the values of calcu- 
lated observables to their values in the full, infinite space [18- 
20, 39, 40], such methods are still in their formative stages, 
especially for the crucial E2 observables. It is therefore partic- 
ularly notable that quantitatively well-developed and robust sig- 
natures of rotation may be observed in the present results. That 
this is possible reflects the distinction between convergence of 
individual observables, taken singly, and convergence of rela- 
tive properties, such as ratios of excitation energies or ratios of 
quadrupole matrix elements. It is these latter relative proper- 
ties which are essential to identifying rotational dynamics and 
which are found to be sufficiently converged to yield stable ro- 
tational patterns at currently achievable N mRX truncations, as il- 
lustrated in Fig. 6 for the K = 3/2 ground-state band of 9 Be. 

From the results of Sec. 3, it is seen that rotational struc- 
ture is pervasive in ab initio NCCI calculations of light nuclei, 
occurring in the yrast and near-yrast regions of all the spectra 
considered for the Be isotopic chain. With suitable extrapola- 
tion methods in place, a salient test of ab initio calculations and 
interactions will then be quantitative prediction of collective ro- 
tational parameters, such as the intrinsic quadrupole moment, 
for direct comparison with experiment. One may observe that 
the present discussion represents a phenomenological rotational 
analysis, in the traditional experimental sense, but of a large set 
of observables taken from ab initio calculations of nuclei. Hav- 
ing full access to the calculated wavefunctions, we may also 
hope to extract information on the collective structure of the 
nuclear eigenstates from other measures of the wave function 
correlations, such as density distributions [20] and symmetry 
decompositions [41]. Natural questions include the origin of 
rotation in the Be isotopes — for instance, the extent to which 
it might arise from relative motion of alpha clusters or perhaps 
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